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A miscut (vicinal) crystal surface can be regarded as an array of meandering but non-crossing steps. 
Interactions between the steps are shown to induce a faceting transition of the surface between a 
homogeneous Luttinger liquid state and a low-temperature regime consisting of local step clusters 
in coexistence with ideal facets. This morphological transition is governed by a hitherto neglected 
critical line of the well-known Calogero-Sutherland model. Its exact solution yields expressions for 
measurable quantities that compare favorably with recent experiments on Si surfaces. 

PACS numbers: 68.35Rh, 5.30Fk 
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Miscutting a crystal at a small angle with respect to 
one of its symmetry planes produces a vicinal surface [Q . 
It often consists of a regular array of terraces separated 
by monoatomic steps. The steps meander by thermal ac- 
tivation but they do not cross or terminate; their density 
is determined by the miscut angle. This picture, the well- 
known terrace-step-kink model |Q , neglects the formation 
of islands, voids, and overhangs on the surface, and is 
hence expected to be valid below the roughening tran- 
sition of the ideal facet. In the simplest approximation, 
such steps are modeled as the world lines of free fermions 
moving in one spatial dimension r and imaginary time 
t, thus taking into account the no-crossing constraint 
through the Pauli principle |^,^. While the free fermion 
model is sometimes a qualitatively satisfactory approxi- 
mation 1^ , it has become clear that interactions between 
the steps can induce phase transitions that change the 
surface morphology drastically ||^. From a theoretical 
point of view, models of interacting fermions are impor- 
tant realizations of two-dimensional euclidean field the- 
ories, some of which are exactly solvable. For example, 
steps with short-ranged interactions can be mapped onto 
the Thirring model or, in the more complex case of re- 
constructed surfaces, onto the Hubbard model 1^. 

Interactions between steps are produced by a variety 
of physical mechanisms 0. For example, elastic forces 
lead to a long-ranged mutual repulsion that decays as 



with the step separation 



Short-ranged in- 



teractions (including all forces that decay faster than 
r~^) can be of either sign. Using scanning tunneling mi- 
croscopy on Cu surfaces, Frohn et al. [|j have found evi- 
dence for step-step attractions that decay over distances 
of a few lattice spacings. In a beautiful series of X-ray 
scattering experiments. Song and Mochrie | [lO| have re- 
cently discovered an important manifestation of attrac- 
tive forces on miscut Si(113) surfaces in equilibrium. At 
sufBciently high temperatures (T>,1300K), a surface of 
miscut angle 9o is a homogeneous ensemble of fermionic 
steps whose local density p(r, t) has the expectation value 
(p(r, t)) — tan 9q = po and somewhat smaller fluctu- 



ations {p{r,t)p{r' ,t')) than expected for free fermions. 
As the temperature is lowered, however, the fluctuations 
increase substantially, until a faceting transition occurs 
at a temperature T*{po) w 1200K. Below that tem- 
perature, the attractive forces cause the steps to clus- 
ter locally. Hence the surface splits up into domains of 
an increased and temperature- dependent step density 
{p{r,t)) = p{T) > po alternating with step-free (113)- 
facets {{p{r,t)) = 0). A critical temperature Tc — 1223K 
is identified from the extrapolation p(Tc) — 0. 

This Letter is devoted to a theoretical analysis of 
stepped surfaces with both long- and short-ranged inter- 
actions. By mapping the step ensemble onto an exactly 
solvable fermion model, it is shown that a long-ranged 
repulsion and a short-ranged attraction can conspire to 
produce a quite complex temperature dependence of the 
surface morphology, including a faceting transition in 
quantitative agreement with the experimental findings 
of ref. jl^. I obtain three temperature regimes charac- 
terized by qualitatively different step configurations (see 
Fig. 1): 






FIG. 1. Typical configurations of non-crossing (fermionic) 
steps coupled by inverse-square and short-ranged forces, 
(a) High-temperature regime [T > T{po)). (b) Critical regime 
{f{po) >T> T*(po)). (c) Faceted regime (T < T*{po)). 

(a) Well above the critical temperature Tc, the steps 
are dominated by the no-crossing constraint and the long- 
ranged repulsion. Hence, they are well separated from 
each other with relatively small fluctuations. Below a 
crossover temperature T{po), the short-ranged attraction 
becomes important, (b) In the critical regime close to T^, 
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the probability of a step being close to one of its neighbors 
is substancially enhanced. This goes along with increased 
step fluctuations and a broader distribution of terrace 
widths, (c) Below the faceting temperature T*{po) < Tc, 
the steps form local bundles of density {p{r, t)) > po. On 
average, the distance between two neighboring bundles 
is larger than the width of an individual bundle. The 
fluctuations of these "composite" steps are smaller than 
those of individual steps. 

Specifically, I consider a system of p fcrmionic lines 
ri{t) governed by the effective action 



S = 



dt 



1 ^ 
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where = dr.i/dt and rij = Ti — rj. The action con- 
tains kinetic terms with a line tension normalized to 1, 
"contact" forces Sa{r) of microscopic range a and 
an "equal-time" approximation ujair) — r^^ to the elas- 
tic interactions for |r| > a The universal properties 
can be expressed in terms of the rescaled coupling con- 
stants go = g/T"^ and /iq = h/T"^. In the limit poa <C 1 of 
small miscut angles, this system can be mapped onto the 
Calogero-Sutherland model |Q, an exactly solvable con- 
tinuum theory well known in the context of the fractional 
quantum Hall effect and random matrix theory. Its two 
branches of solutions are labeled by the parameter 



A±(5o) 



1 ± VI + 4.go 
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The strong temperature dependence of the surface mor- 
phology described above is shown to arise from crossover 
phenomena between these two branches of solutions. In 
the high-temperature regime, the surface is governed by 
the solution X^{g/T'^). At the critical temperature Tc 
(implicitly given by h/T^ = (a, g/T^), where (a, go) 
is a nonunivcrsal function), the surface scales according 
to the solution A^((jfg) (as long as .gg = g/T^ < 3/4; be- 
yond that point, this branch of solutions ceases to exist). 
The solution A~ (.g§) also determines the singular density 
dependence of the crossover temperature T{po) and the 
faceting temperature T*{po) (i.e., of the size of the criti- 
cal regime). Hence, a number of observables at different 
temperatures are predicted in terms of the single nonuni- 
vcrsal parameter g^, which allows direct comparison with 
experiments (see the discussion at the end of this Letter) . 

The solutions of the Calogero-Sutherland model la- 
beled by < A < cxD are known to describe a line of 
Luttinger liquid critical points ||l^ that contains the self- 
dual point A = 2, the free fermion point A = 1 and 
the Kosterlitz-Thouless point A = 1/2. Notice, however, 
that for repulsive long-ranged forces {g > 0), the root 
A (.go) is negative. Solutions of the Calogero-Sutherland 
model with A < have not been discussed before as 
they were deemed unphysical. The solutions labeled by 



0>A>— 1/2 form a new line of Luttinger liquid critical 
points; faceting on vicinal surfaces seems to be their first 
realization. This line is the analytic continuation of the 
line < A < oo beyond the free boson point A = (see 
Fig. 2), and it terminates at its Kosterlitz-Thouless point 
A = —1/2. There is another closely related physical man- 
ifestation of the solutions X~(go) in the particular case 
p = 2, where the action ([T]) is a model for two interfaces 
in a two-dimensional system in the so-called intermediate 
fluctuation regime p^ . In these systems, the well-known 
line of wetting critical points p^ , p^ turns out to corre- 
spond to that branch of solutions. The wetting transi- 
tion is of second order for A > —1/2, but of first order 
for A < -1/2 lll,|ll. 
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FIG. 2. The two branches X'^{go) and \~ (go) of the 
Calogero-Sutherland model. On vicinal surfaces, they gov- 
ern the high-temperature regime and the faceting transition, 
respectively. The solutions < A < cxd and > A > — 1/2 
form two distinct lines of Luttinger liquid critical points. Spe- 
cial points are A = 1 (free fermions), A = (nonrelativistic 
free bosons) and A = ±1/2 (Kosterlitz-Thouless points). 

To derive these results, it is convenient to regard the 
ensemble of steps as a many-body quantum system. The 
Hamiltonian of this system. 



H 



dri;\r,t)d^i;{r,t)+gona{t) + ho^ait), (3) 



acts on p-particle states; its form is determined by 
the action (^. ip and are anticommuting fields; 
nait) = Jdrdr'p(r,t)u;a{r - r')p{r',t) and $^(0 = 
J drdr' p{r,t)Sa{r — r')p{r' ,t) are the long- and short- 
ranged interactions written in terms of the density op- 
erator p{r, t) = ip''{r, t)ip{r, t). In a system of finite width 
L with periodic boundary conditions and the periodic po- 
tential uja{r) = (tt^/X^) sin^^(7rr/i) (for a < r < L — a), 
the two-particle ground state takes the exact form 



^'2(ri2) = sin^(7rri2/L) 



(4) 



(for a < ri2 < L — a), with A given by Eq. (^). A contact 
potential of fixed strength ho = h^{a,go) is required to 
match the wave function (Q) with the fermionic boundary 
condition 5-2 (0) = *2(i) = ||l||. In the limit a 0, 
the p-particle ground state is the simple product 
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(5) 



provided A > —1/2 (otherwise this wave function is not 
normaUzable). An integrable continuum model emerges, 
known as the Calogero-Sutherland model [Q. 

Hence, for a given value —1/4 < go < 3/4, the Hamil- 
tonian defines two different continuum theories, corre- 
sponding to the branches X^{go) in Eq. (H). The stability 
of these solutions in the thermodynamic parameter space 
(50 J ^0) can be studied perturbatively, using the methods 
of refs. |p^|-pl[| (where the reader is referred to for more 
details). The expansion of the dimensionless p-particle 
ground state energy Fp = L'^Ep about the branch point 
A* = 1/2 has the form 



Fp{g'o,K)~Fp{0,Q) 



(-1) 



M+N 



M\N\ 



M+N>1 



,M,,N 
i'M,N 90 "0 



with 5o = 50 + 1/4, h'^ = ho - h„ (a, -1/4), and 

, M+N M M+N 
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M,N = Po ' 



j n dt.(ni^a(i.) n *"(^«)) 

i=2 1=1 i=M+l 



otherwise. The brackets (. . .) denote connected expec- 
tation values in the unperturbed p-particle ground state 
(^ at A = 1/2. In the limit a ^ 0, the coefficients Fm,n 
develop logarithmic singularities. With an appropriate 
normalization of the operators $a and f2a, one finds at 
the lowest orders 

Fo^2 = Po^ J dt ($a(0)$a(i)) = ^SPo^i'^a) + 0(s°) , 

and hence 

Fp{gi h'„) - Fp(0,0) - ~po'{^a){h'o + sig'o - h'^')) 
+0{s°,g'Q^,g'„hQ,hQ^) , 

where s = — log(poa)- Up to this order, the singularities 
can be absorbed into the renormalized coupling constant 
hj^ = ft,Q(l — s/i/f) + sgQ + . . ., while no renormalization is 
needed for go- This leads to the parabolic flow equations 

d d 1 2 

^ go = , —Hr = ^ + .go - /ifi , (6) 

which are independent of p. They have first been ob- 
tained by functional renormalization group methods for 
p = 2 1 17 1 and have been derived for arbitrary p in 
refs. ||2^,23|. It is possible to check that Fi^ and ^0,2 
contain the only primitive singularities of the perturba- 
tion series. The renormalization group equations are 
thus exact to all orders in a minimal subtraction scheme. 

There are two lines of fixed points, — ±-y/l + 4(7o/2. 
The renormalization group eigenvalue of temperature 
variations (or variations of hu), 



1 d 



2 dhR 



fdhi 



V ds 



VI + 450 



(7) 



also governs the scaling of the contact operator e.g., 
{^a)'^igo) ~ Po^'^ 'y3»)) ^ xhis is precisely the scaling 
of (^'p|<i>Q|\E'p) obtained from the exact solution (||) with 
A given by (||). The two lines of fixed points ^.^(50) 
can thus be identified with the two branches of solutions 
\^{go) of the Calogero-Sutherland model. 

At the fixed points A~((7o), temperature variations are 
a relevant perturbation. (For p = 2, these fixed points 
govern the wetting transition at the critical temperature 
Tc [|3-) Above Tc, Eq. (|) yields a crossover to the 
stable branch of solutions \^{go)- The crossover temper- 
ature is given by f{po) - - pl^'^'^^'K Below T^, the 
renormalized coupling tends to —00 under the fiow 
(^), indicating an instability of the step ensemble with 
respect to the formation of local bundles. These bundles 
have a characteristic line density 



p{T) ~ (Te-T)i/2!/"(ffS) 



(8) 



The clustering becomes visible if p{T) > po, i-e., for T < 

T*{po) with Tc-T*{po) ~ pY}^°^- On average, a bundle 
consists of n w 20 lines and has a width n/p{T). 

Two neighboring bundles at a typical distance n/ po have 
an exponentially small overlap ~ exp(— p(T)/po); they 
can thus be approximated as stable composite steps with 
nearly step-free facets in between. 

From the above discussion, it is clear that the integra- 
bility of this system is tied to scale invariance at distances 
a <^ r <^ Pq^ . Along the crossover between the critical 
and the high-temperature regimes, the two-particle wave 
function does not have the simple power-law asymptotics 
^'2(''') ~ l?'!'^ as in (^, and consequently, the product 
ansatz (|^) breaks down. What happens to scale invari- 
ance at distances r 3> Pq ^? It has been shown that any 
solution of the Calogero-Sutherland model with A > de- 
scribes a Luttinger liquid: it belongs to the universality 
class of the Gaussian model with action 



= ^ / (VMr))^d2r 



(9) 



and stiffness 7 = A/2 |T5|. Here r = {r.vpt), where 
vp = 2ttXpo is the Fermi velocity, and h(r) is a coarse- 
grained surface height variable. Using the Bethe ansatz, 
one finds the low-lying finite-size excitations AEe^m = 
2TTVFXe,m/ L in terms of the scaling dimensions Xe^m = 
(e^7 -I- rr? Ij)j2 of the Gaussian vertex operators Oe.m 
(e, TO e Z) It is then easy to show that any solution 
with A < is also a Luttinger liquid with 7 = —A/2, 
since the transformation A ^ — A acts as the symmetry 



on the Gaussian operator algebra. Thus 



conformal field theories with central charge c = 1 govern 
the steps at Tc and in the high-temperature regime - and 
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thus along the entire crossover by virtue of the c-theorem 
psf . It follows that the system is a Luttinger liquid at 
all temperatures T > pGf . This property extends 
to the entire critical regime, where 7 can be written in 

scaling form, 7(r,po) = r((r - Tjp^^^"^^"^). Below 
T*{pa), one expects an effective action similar to for 
the composite steps; the system is then still a Luttinger 
liquid. 

We now turn to some experimentally measurable con- 
sequences of this theory. 

(a) The ground state wave function yields immedi- 
ately the step density correlation function 

{p{0,t)p{r,t)) ^^l{r)^r'^ (10) 

in the limit |r| ^ Po^ ^ where multi-particle effects can 
be neglected. Two steps at distance r from each other 
enclose a terrace of width r if no further steps are in 
between; this condition is also negligible for r ^ Pq^. 
Hence the terrace width distribution, which can be mea- 
sured by surface scanning techniques, has the same short- 
distance tail ^ r^^ In the high-temperature regime 
(A = A+(g/T^) > 1), short terraces are rare, while at 
(A — X^igo) < 0), they are abundant. 

(b) The surface stiffness 7 is measurable as universal pre- 
factor of the height difference correlation function 

C(r) ^ i((/i(0,t) - h{r,t)f) 7-'log|por| + ... (11) 

on scales |r| 3> Pq^ ■ 7 shows a characteristic temperature 
dependence. Its high-temperature asymptotic value is 
7 = 1/2. As the temperature is lowered, it first increases 
as 7 = X+{g/T^)/2 = X+ {g^T^ /T^)/2, then decreases in 
the critical regime (taking the value 7 = — A~(f;g)/2 < 
1/4 at Tc), and again sharply increases below T*{pq) to 
values 7 > 1/2 for composite steps. Most aspects of this 
pattern have been observed jl^, but the measurements 
are not yet conclusive in the high-temperature regime. 
In the critical regime, the data have been fit to a power 
law 7 ~ (T-r,(po))-" with T,(po) < T*{po) (lO), which 
corresponds to a singularity in the scaling function F. 

(c) The temperature dependence of the cluster density 
("^ agrees with the measurement p{T) ^ (7;„r)0-42±o.io 
m ifff§«3/4. 

(d) The temperature differences f{pn) ~ T^, - T*{po), 

and Tc — Ts{po) depend on the step density as p^^ 
This is also consistent with the data for w 3/4. 

In summary, the Calogero-Sutherland model has been 
applied to interacting steps on vicinal surfaces. The ther- 
modynamic complexity of this system arises from the in- 
terplay of the two branches of integrable solutions. It will 
be of interest whether this mechanism also plays a role 
in other realizations of the Calogero-Sutherland model. 
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